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The structure of the neutrino mixing matrix is indicative of an underlying family symmetry that 
interrelates the three generations of fermions in the Standard Model. We systematically scan 
the parameter space of 76 discrete non-Abelian family symmetries and construct all models 
with the Standard Model particle content and up to three flavon fields where we include non- 
renormalizable interactions of mass dimension five and six. We find that of the 76 groups that 
we considered, 44 groups can accommodate models that are consistent with experiment at 3a, 
and 38 groups can have models that are tribimaximal. One immediate consequence is that A4 
is not "special", but should be considered on equal footing with other groups such as T-j that is 
the smallest group for which we find tribimaximal mixing, and Ti^ that has the largest fraction 
of TBM models. We present the details of a model with 612 = 33.9°, 023 = 49.1°, 613 = 5.1° to 
show that a non-zero 613 can easily be accommodated. 
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1. Introduction 

The discovery of neutrino masses is one of the most exciting developments in physics during 
the last decade. Unlike the CKM matrix in the quark sector, the mixing matrix ?7pmns contains large 
angles, and until very recently its form was well-approximated by a pattern that was subsequently 
dubbed tribimaximal [1]. The results from the T2K [2] and MINOS [3] collaborations have called 
this paradigm into question, and global fits to neutrino oscillation data now indicate that On ^ 
with about 3a significance [4]. 

The (near) tribimaximal mixing pattern has given rise to speculation that there might exist a 
symmetry that inteiTclates the three generations of leptons in the Standard Model. The transforma- 
tion properties of the leptons under this so-called family symmetry then determine the Lagrangian 
and ultimately the mixing pattern. For a recent review on finite symmetries and their application 
to neutrino physics, see Ref. [5]. In the past, models based on S3, D4, A4, S4, A5, T', A(27), to 
name a few, have been particularly popular, and it has been ai^gued that A4 is inherently connected 
to tribimaximal mixing (TBM). Here we describe a systematic approach to the problem of finding 
the family symmetry that may best describe the observed pattern of mixing angles [6]. To that end, 
we use GAP [7] for the representation theory of discrete groups and an algorithm due to van den 
Broek and J.F. Cornwell [8] to calculate the Clebsch-Gordan coefficients for any finite group. The 
availability of such powerful algebraic tools and the fact that the Lagrangian is fully determined by 
the underlying symmetries of the theory render the calculation of neutrino mixing angles an ideal 
problem for a computational approach. 

2. Algorithm for the Scan 

There are 1,048 discrete groups with less than or equal to 100 elements. Since we want to 
interrelate the three generations of leptons, we will put the lepton doublets L^, L^, into one 
and the same irreducible triplet representation of the family symmetry. This leaves us with 90 
groups that have three-dimensional irreducible representations. It turns out that for 14 groups a 
comprehensive scan would take too long, so in the end we are left with 76 groups [6]. For the right- 
handed leptons e, IJ.,T and the two Higgs doublets /j„ and h^i, we assume that they transform in one- 
dimensional (not necessarily singlet) representations. The choice for the Higgs fields is motivated 
by supersymmetry, whereas the assumptions for the right-handed leptons should eventually be 
generalized. We spontaneously break the family symmetry by giving vacuum expectation values 
(vevs) to flavon fields (pj, (ps, ^ that are gauge singlets, but may transform in any representation of 
the family symmetry. We include in the Lagrangian all terms of mass dimension less than or equal 
to 6. In scanning over the vevs for the flavon fields we have not been fully general, but restricted 
ourselves to those with entries or 1 (times a constant). We agree that we are sampling only a 
discrete part of the parameter space, but choosing continuous values is computationally not viable. 
Moreover, that the vevs are all of the same order of magnitude is not an unreasonable assumption. 

3. Discussion of Results 

In Fig. 1 , we summarize the results of our scan. Not all of the groups that we have considered 
have names in the conventional sense, so in the following we will denote them by their so-called 
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GAP IDS [6]. E.g. 0(12,3) is the alternating group A4, i.e. it is the third in the list of all groups of 12 
elements. Out of the 76 groups that we have scanned, 9 (12%) have only singular mass matrices. 44 
groups (58%) lie in the 3a interval^, and 38 (50%) are even tribimaximal (for at least one vacuum 
configuration, respectively). We find that the smallest group that can realize TBM is Tj, and the one 
with the largest fraction of TBM models is T^. Regai^ding whether there is a special connection 
between A4 and TBM, we find that taking the fraction of TBM models as a criterion, T] and Tn 
are as promising, if not more, than A4. Furthermore, only 16 of the 35 groups containing A4 as a 
subgroup (red and green labels on the jc-axis) realize TBM, while there are 22 groups which do not 
have A4 as a subgroup but do realize TBM. 
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Figure 1: The number of models per symmetry group. The green and the red bars indicate the number of 
models that for given vevs lie within the 3(7 interval or are tribimaximal mixing, respectively. 

In Fig. 2, we show the distribution of the mixing angles for one particular group, namely 
Tij. In the first histogram we count the number of vacua that realize du as given on the x-axis, 
where 623 and 613 are restricted to their 3a intervals. Analogous statements hold for the two other 
histograms. The maxima in all three histograms correspond to tribimaximal mixing, and we obtain 
the almost unique prediction 023 = 45°. The most preferred value for 613 is 0°, but non-zero values 
are also allowed in the 3 a band. 
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Figure 2: Number of vacua with family symmetry C5(39, 1) — Ti^, that give the mixing angles on the .x-axis. 



4. A Model with 613 ^ 

We present a specific model [6] where 613 / and all mixing angles are compatible with 
experiment [4] at 2a. For the family symmetry, we take 6(36,11) = A4 x Z3 and assign the 
following family charges to the Standard Model particles and flavon fields: 

{L,e,^,z,hu,hd,(pT,(P s,^) ~ (3, 1, 1'", 1(4), 1', l(^), l(5\3",3"). (4.1) 

'Note that in this section we refer to the experimental limits that we used in our original work [6]. 
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The Lagrangian is given by all the terms that are invariant under both the gauge and the family 
symmetries and are at most of mass dimension 6: 

ct LLhuKcpr + csLLhuhu(ps+c^ LLhuK^ + d^^^ Leh^cps + d^^^ Leh^^ 

Here, we have only listed those terms that will contribute to the mass matrix. The c's and J's 
are some coefficients that can be related to the Yukawa couplings in a basis where the charged 
lepton mass matrix is diagonal. After contracting the family and SU(2) gauge indices and giving 
vevs to the flavon fields along the directions {(pj) = vj, {(ps) = ivs,^,vs), {^) = (v^,v^,v^), we 
obtain the mixing angles du = 33.9°, 623 = 49.1° and 613 = 5.1° that all lie in the 2a interval 
of the experimentally determined values. We have assumed that all vevs are of the same order of 
magnitude (the absolute value of this scale does not enter the calculation of the mixing angles). 

5. Conclusions & Outlook 

In this work, we pursued a systematic search for models of neutrino flavor that reproduce the 
observed mixing pattern. We find that half of the 76 flavor groups that we considered can give rise 
to tribimaximal mixing. As a consequence, we cannot affirm that A4 is "special", in particular as 
compared to Ty or Ti^. Following the recent hints of a non- vanishing 613 we presented a model 
where all mixing angles are compatible with experiment at 2a. In a next step, we will extend our 
analysis to include neutrino masses and the stabilization of the flavon vacuum expectation values. 
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